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We report on the temperature dependence of the mobility, fi, of the two-dimensional electron gas 
in a variable density AlGaN/GaN field effect transistor, with carrier densities ranging from 0.4 x 10 12 
cm~ 2 to 3.0xl0 12 cm -2 and a peak mobility of 80,000 cm 2 /Vs. Between 20 K and 50 K we observe a 
linear dependence p^c = <*T indicating that acoustic phonon scattering dominates the temperature 
dependence of the mobility, with a being a monotonically increasing function of decreasing 2D 
electron density. This behavior is contrary to predictions of scattering in a degenerate electron 
gas, but consistent with calculations which account for thermal broadening and the temperature 
dependence of the electron screening. Our data imply a deformation potential D = 12-15 eV. 
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In high quality AlGaN/GaN heterostructures the 
temperature dependence of the mobility of the two- 
dimensional (2D) electrons is determined solely by 
acoustic phonon scattering over a wide range of low 
temperatures 0, 0- In this regime the inverse electron 
mobility, which is proportional to the scattering rates, 
varies linearly with temperature with the slope depend- 
ing only upon the 2D electron density, and the values 
of the deformation potential and piezoelectric constants 
which characterize the strength of the electron-phonon 
interaction. Therefore, by measuring the temperature 
dependence of the mobility across a range of densities, 
and using literature values of the piezoelectric constants, 
both the poorly known deformation potential and the 
density dependence of electron scattering rates due to 
acoustic phonons may be determined. 

The density dependence of acoustic phonon scatter- 
ing in GaN heterostructures has not been extensively ex- 
plored. To date, the only data on acoustic phonon scat- 
tering in the AlGaN/GaN 2D electron system (2DES) 
are from a small collection of fixed, moderate density 
samples 1] . However with the recent development of high 
quality, high mobility AlGaN/GaN field effect transistors 
(FETs), it is now possible to study transport phenomena 
across a wide range of densities in a single device By 
applying a voltage to the gate of the FET, the 2D elec- 
tron density, n2_o, may be continuously varied with the 
advantage that all other system parameters (e.g. the dis- 
tribution of fixed imperfections) remain constant. We 
have utilized such a high mobility FET, together with 
a set of fixed density heterostructures, to systematically 
explore the density dependence of acoustic phonon scat- 



tering rates over more than an order of magnitude in den- 
sity. This range is much broader than in previous work 
and so allows an effective comparison with theories of 
electron scattering by deformation potential and piezo- 
electric phonons. We find that existing calculations of 
acoustic phonon scattering in degenerate GaN 2DESs are 
insufficient to explain our data quantitatively. However, 
the main features and trends can be reproduced if the 
temperature dependence of electron energy and screen- 
ing are taken into account. 

In this study we have utilized a FET along with a set 
of fixed density heterostructures. The FET device has a 
peak mobility of 80,000 cm 2 /Vs, and a continuously tun- 
able n-2D that ranges over 0.4 — 3.0 x 10 12 cm -2 , among 
the lowest densities measured in GaN. The fixed density 
samples have = 1 — 7 x 10 12 cm -2 , and low-T mobili- 
ties of 8, 000 — 30, 000 cm 2 /Vs. All devices were grown by 
molecular beam epitaxy (MBE) on GaN templates pre- 
pared by hydride vapor phase epitaxy (HVPE). In the 
FET, the HVPE GaN is 40 ^m thick and is grown atop 
a sapphire substrate. The HVPE GaN is known to have 
a low density of threading dislocations of approximately 
10 8 cm -2 , and has been compensated with Zn at 10 17 
cm~ 3 to suppress any residual conductivity |3j. A layer of 
1.5 fim nominally undoped MBE GaN is grown on top of 
the HVPE GaN, and followed by 160 A of Al . 06 Ga .94N. 
The capping layer consists of 30 A of MBE grown GaN. 
The fixed density samples were grown similarly, but with 
different AlGaN and GaN thicknesses in order to vary 
112D- To make the FET, a 0.5 mm long and 80 /im wide 
Hall bar was patterned by photolithography, followed by 
a 100 nm deep etch in a Cl-based plasma. Three voltage 
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FIG. 1: Typical experimental mobility vs. temperature, 
showing regions dominated by different types of scattering. 
Inset: inverse mobilities plotted against temperature for two 
different densities in the FET; solid lines are a guide to the 
eye. 

contacts are located on either side of the Hall bar, and an 
insulated gate is created on the FET surface by deposit- 
ing 10 nm Ni and 100 nm Au on top of 100 nm of Si02- 
The fixed density samples consist of 4x4 mm squares 
with 8 contacts placed symmetrically around the edge. 
The contacts in all structures consist of a Ti/Al/Ni/Au 
metal stack annealed at 800° C for 30 s. 

We performed low field electron density and mobility 
measurements in the FET and fixed density samples in 
the standard Hall and Van der Pauw configurations. Fig. 
1 shows a typical trace of mobility vs. temperature, in 
which regions limited by different scattering mechanisms 
are marked. Mobility data were taken between 300 and 
1.7 K in most devices. Between 20 K and 50 K, a clear 
linear relationship is observed when the data are plotted 
as inverse mobility vs. temperature; see the inset of Fig. 
1. In some traces (e.g. the top one in the inset) the data 
depart from linearity as T — > 50 K; we will return to this 
observation below. In common with previous authors 1, 
Q we perform a linear fit to the inverse mobility using 



1 

M 



1 

Mo 



1 

f^ac 



l 

Mo 



aT 



(1) 



where a is the slope of the fit (a = d(l/ii)/dT ), [i ac is the 
component of the mobility due to acoustic phonon scat- 
tering, and the extrapolated zero-temperature mobility 
Ho accounts for all T-independent scattering mechanisms 
such as from remote impurities or charged dislocations. 
Fig. 2 shows the slope a as a function of ii2D- The sig- 
nificance of a lies in its proportionality to the acoustic 
phonon scattering rate, l/r ac = (e/m*) aT. Therefore, a 
determines the density dependence of the electron scat- 
tering rates due solely to acoustic phonons. 
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FIG. 2: a vs. 2D electron density. Data shown is from the 
FET (■), fixed density samples (♦), and Ref[l] (A). Calcu- 
lated fits are for two values of the deformation potential (D = 

12 and 15 eV): — for degenerate 2DES; for non-degenerate 

2DES, with ridepi = 0; and ••• for non-degenerate 2DES, with 
ndepi = 6 x 10 11 cm" 2 . 



The upward trend of a at low n.2D in Fig. 2 is striking. 
It departs strongly from predictions based on a treatment 
of acoustic phonon scattering in a screened, degenerate 
2DES (solid lines) 0. We see a increase monotoni- 
cally by nearly a factor of two over an order of magni- 
tude decrease in density for ii2D < 2 x 10 12 cm" 2 , and 
show no particular density dependence for n2D > 2 x 10 12 
cm" 2 . Data marked as A are due to Knap et al 1], which 
agree with our high density data, n2£> ~ 5.5 x 10 12 cm" 2 . 
However the lower density data of Ref [1] fall considerably 
below ours although they are close to the degenerate the- 
ory. The source of this discrepancy is unknown, but is 
unlikely a result of the differences in scattering from fixed 
sources. Such an effect is included in the term /ioj which 
for n 2D = 2- 3x 10 12 cm" 2 is 70,000 cm 2 /Vs in our 
FET and « 20, 000 cm 2 /Vs in the fixed density samples; 
the lower density samples of Ref[l] have Ho « 60, 000 
cm 2 /Vs. 

A qualitatively similar trend in the density dependence 
of a has been previously observed in the AlGaAs/GaAs 
2DES by Harris et aZ@. Their data also conflicted with 
predictions of phonon scattering in a screened, degen- 
erate 2DES. This was resolved by Kawamura and Das 
Sarma 7j, who noted that at lower n.2D the Fermi tem- 
perature, which is proportional to n2D, becomes compa- 
rable to the highest temperatures in the experiment (e.g. 
Tp = 41 K for ii 2 _d = 1 X 10 11 cm" 2 in GaAs; Harris 
et al took data up to 40 K). The 2DES now ceases to 
be degenerate and the simple linear temperature depen- 
dence of is modified. To account for this, Ref [7] cal- 
culates acoustic phonon scattering in a non-degenerate 



3 



2DES by including thermal broadening of the electron 
distribution function and the T dependence of electron 
screening. These additional T dependencies tend to de- 
crease the mobility as T — > Tp, resulting in a nonlinear 
increase of as T increases. Therefore, a linear fit to 
fi^ 1 using eq. [JJwill create a slope (a) greater than if non- 
degeneracy effects were not present. This is clearly visible 
in the inset to Fig. 1: the top trace, for 112D = 0.66 x 10 12 
cm~ 2 (Tp = 80 K), departs from linearity toward 50 K; 
whereas the bottom trace, with n 2 £> = 1.67 x 10 12 cm~ 2 
(Tp = 202 K) is linear over 20 - 50 K. 

Guided by the similarity between the GaN and GaAs 
data, we adapted the calculations by Priced and Kawa- 
mura et al 7] for the GaN system, using the values of 
the deformation potential, D, and the depletion charge 
density, n^epj, as adjustable parameters. All scattering is 
with acoustic phonons due to the deformation potential 
and the piezoelectric field. In calculating the contribu- 
tion from the latter, we used the form factors appropri- 
ate for scattering of 2D electrons by bulk phonons in 
a wurtzite crystal @. We used Q = 3.7 x 10 11 N/m 2 
and C t = 1.3 x 10 11 N/m 2 for the longitudinal and 
transverse elastic constants^, with /i 33 = 0.65 C/m 2 , 
^31 = —0.37 C/m 2 and his = —0.33 C/m 2 for the piezo- 
electric constants |sL TToj . 

In Fig. 2 we have plotted the results of our calculations 
for comparison with our data. The solid lines show the 
theory for the degenerate case with deformation potential 
D = 12 and 15 eV, and zero n.d ep i- The dashed and dot- 
ted lines are calculated using the non-degenerate theory 
of Kawamura et al for D = 12 and 15 eV, and various 
v-depi- Clearly these reflect the data much better than 
the degenerate theory, in particular by qualitatively cap- 
turing the rise in a at lower densities. The best fits are 
achieved by using non-zero values for Hd ep i, although no 
single set of values for D and iid ep i reproduces a over the 
measured range of n2£>. We achieve the best results at 
lower densities for nd ep i = 6 x 10 11 cm~ 2 and D = 15 eV 
(upper dotted line), while for n2£> > 2 x 10 12 cm" 2 , the 
better fit is with D = 12 eV. These values of D are con- 
siderably higher than the D = 8.3 eV calculated for bulk 
wurtzite GaN^lJ, analogous to the situation in GaAs for 
which D = 7 eV in bulk, and 13 eV in 2D scattering^. 



In comparing our data with our calculations we note 
that piezoelectric phonon scattering becomes important 
at low densities [jj, but the piezoelectric constants h x of 
wurtzite GaN are not well known. Prior to the work of 
Kawamura et al, Harris et al noted that in GaAs the rise 
in a at low densities could be reproduced by enhancing 
the value of hi4 (for cubic GaAs) by a factor of 1.50. 
Thus the rise of a in the GaN 2DES may be partly due 
to larger values for the h x than reported in the litera- 
ture. Also, the value of rid ep i required in the best fit 
seems rather high. Assuming pinning of the Fermi level 
at mid-gap in the Zn-compensated HVPE GaN, the elec- 
tric field across the undoped MBE GaN is equivalent to 
only Udepi = 1.2 x 10 11 cm -2 . Finally we note that at 
high n2_o, the Bloch-Gruneisen (B-G) effect will reduce 
phonon scattering for T < Tb_g0I- We have repeated 
our analysis on fi ac vs. T for — 10 13 cm~ 2 from 
Fig. 1 of Ref[12], for which T B _ G = 94 K, and find 

that a is increased by w 9% in the B-G regime. With 

1/2 

Tb-g ~ ^2p> 1 thi s effect rapidly diminishes with decreas- 
ing density and is not involved in the low density behavior 
of a. 

In conclusion, we have measured the temperature de- 
pendence of the 2D mobility over a large range of electron 
densities in high quality AlGaN/GaN heterostructures. 
We compared these results with calculations of electron 
scattering by acoustic phonons, and found a qualitatively 
good fit only by accounting for additional temperature 
dependencies of the electron gas that arise when the sys- 
tem becomes non-degenerate at low densities. We de- 
rive a value for the GaN acoustic deformation potential 
between 12 and 15 eV, with the fit improving with an 
ad hoc, non-zero depletion density. It remains unclear 
whether the present theory, even when adjusted for more 
precise future piezoelectric constants, will be able to ac- 
count for the experimental data. 
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